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Abstract 

Blisters arc commonly observed when an adhesive sheet is carelessly deposited on 
a plate. Although such blisters are usually not desired for practical applications, we 
show through model experiments on angular blisters how material properties can be 
deduced from height profile measurements. In particular the typical curvature of the 
crests is found proportional to an elasto- capillary length that compares the bending 
stiffness of the sheet with adhesion energy. In addition, the radius of the tip allows to 
estimate the product of this length with the thickness of the sheet. The relevance of 
these results to realistic random configurations is finally confirmed. 

1 Introduction 

Placing a sticker on a flat surface is usually a delicate task since blisters form very 
easily. Air may indeed remain trapped between the adhesive sheet and the substrate 
leading to bubble-shaped blisters, while linear folds and wrinkled structures may appear 
when mismatching adhesion fronts are brought together. The consequences of blisters are 
important in a wide range of industrial applications involving the deposition of an adhesive 
sheet on a panel where such defects are crucially undesired. Similar blisters are also 
observed at nanoscale when graphene sheets are deposited on substrates (XJ [2J El S] - Beyond 
deposition processes, blisters are also commonly observed when materials covered with thin 
films delaminate. Indeed differential thermal expansion, differential swelling or residual 
stresses due the coating process can lead to the formation of wrinkles [5] and eventually 
to the debonding of the film if these stresses are compressive [6j |8j [TU] . While most 
studies have focused on predicting the onset of buckling with the aim of avoiding disastrous 
defects in thin film coatings, recent experiments have shown how surfaces can be micro- 
patterned with controlled wrinkles and blisters |5j HT] , with potential applications spanning 
from optics [12J to stretchable electronics |13[ [HJ US]- Another interesting application of 
wrinkles and blisters, following studies on the blister test |16|I17|. is to use their geometrical 
characteristics to probe the mechanical properties of thin sheets as well as the energy of 
adhesion [18, 19, 20j. Motivated by these novel metrology techniques, we propose to study 
which relevant information may be deduced from monitoring the shapes of the blisters 
obtained when a thin adhesive sheet is randomly deposited on a flat substrate. We first 
describe a model experiment complementing, in a different regime, the work from Chopin et 
al. |21j . In this experiment adhesion is due to a film of wetting liquid preliminary coated on 
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the substrate and a single angular blister is produced by compressing one edge of the sheet. 
Results from numerical simulations are in addition compared with experimental data. In a 
second section, the outputs of the simplified configuration are applied to complex blisters 
obtained by depositing carelessly a thin sheet on an adhesive (non slippery) substrate. 
The different physical parameters deduced from monitoring the curvature of the crest and 
the radius of the tip of the blisters are finally compared with their estimates from other 
techniques. 

2 Model experiment with a single blister 

Consider an adhesive film carelessly thrown over a substrate. The parts of the falling 
adhesive sheet that first touch the substrate immediately stick to it. Around these contact 
regions, adhesion zones develop and are generally not compatible with the other patches, 
which leads to the formation of blisters. We usually observe that the blisters often lie along 
a straight direction, with a slowly varying width and a finite length L. In this section we 
study a model experiment on isolated angular blisters with such properties. 

2.1 Experimental setup and typical observations 

Thin sheets of biaxially oriented polypropylene (Innovia films) are laid down a rigid 
plate covered with a layer of ethanol. The Young's modulus and Poisson's ratio of the 
polymer are E = 2.6 ± 0.2 GPa 1 []and v = 0.4, respectively. Three film thicknesses were 
selected for our experiments: /i=15, 30 and 50 fim. Due to the surface tension of the liquid 
(7 = 22.4 mN.m -1 ), the sheets adhere to the plate (sheets and plate are totally wetted by 
ethanol). As sketched in figure [l^, a controlled transverse displacement is imposed at one 
edge, leading to the formation of a single blister that propagates through the sheet up to a 
finite distance L from the compressed edge. The shape of the observed structures evolves 
slightly with the thickness of the liquid layer, but a quasi-steady regime is achieved as the 
liquid progressively evaporates. Our measurements are conducted in this regime, where 
the shape is independent of the remaining volume of liquid. 



Figure 1: (a) Experimental setup: a thin sheet is deposited over a rigid plate coated with 
ethanol. A transverse displacement AIq is imposed at one edge with a caliper symbolized 
by the arrows, leading to the formation of a single blister. The shape of the blister is 
defined by its local height d and width A associated with a local transverse displacement 
Al. (b) A laser sheet tilted by an angle 9 is used to measure the height profile of the 
blister. 

1 The range given for Young's modulus corresponds to the anisotropy due to the process of fabrication. 
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The planar deformation of the sheet is measured using a standard image correlation 
technique (DaVis software from LaVision): the thin films are sputtered with black paint 
and a picture taken after displacement is compared with a reference image of the non- 
deformed film. The cross-correlation of the two pictures leads to the planar displacement 
field of the sheet (arrows in figure |2^). Comparing the motion of the sheet on both sides 
of the blister gives a measurement of the transverse displacement Al along the blister as a 
function of the distance to the edge X e d ge . We observe a linear variation of this transverse 
displacement with X e d ge (figure [2^) , which shows that the opposite sides of the thin sheet 
separated by the blister undergo a relative rotation motion (see also arrows in figure pk) 




around the tip of the blister. Indeed, both patches move as rigid bodies (a composition of a 
translation and a rotation) with respect to each other. At the tip of the blister, this relative 
motion is zero (Al = 0), so that both patches rotate around it. Since the localization of 
the center of rotation is very accurate from these measurements, we define the length L 
of the blister as the distance of this center to the edge. In the following of the paper we 
choose this center as the origin of our coordinate system (with x = L — X e d ge ). 

In our experiments, the length L of the blister was not reproducible, and we do not 
give a prediction for this parameter in this paper. Indeed, the propagation of the blister is 
stopped by secondary wrinkles and blisters forming in front of the tip. As a matter of fact, 
we expect the rotational motion of the thin sheet to generate ortho-radial tension ahead 
of the tip, inducing radial compression (due to Poisson's ratio), which in turn leads to 
the formation of wrinkles or secondary blisters perpendicular to the main angular blister. 
These smaller structures are observed in the vicinity of the tip and are more numerous 
for thinner sheets 2 . [] We focus now on the profile of the blister along its axis. The width 
A, the amplitude d and the transverse displacement Al are measured as a function of the 
distance to the tip along the axis x, for a given displacement AIq at the edge. The height 
profile is obtained by scanning a laser sheet tilted by an angle 9 along the blister (see figure 
[lj>). After calibration, the in-plane displacement of the laser sheet seen from above gives 
the out-of-plane displacement |14| . Thus, the shiny line observed at the intersection of the 
laser sheet with the blister directly gives the local profile (figure [2^,) . The amplitude and 
width are measured by fitting the profile by a cosine function z(y) = d/2[l + cos(27ry/A)], 
since this form is classically observed with beams slightly bent. The variations of d(x) and 
A(x) are respectively displayed in figures^ and[2]l. 

Plotting the profiles of different sections of the same blister in stretched coordinates, 
with z and y respectively normalized by x 2//3 and x 1//3 , reveals a self-similar shape (figure^. 
However, contrary to the case of the liquid blister test studied by Chopin et al. |21| . the 
profile is not conical (which would correspond to both A and d being proportional to the 
distance to the tip x). The aim of the following section is to rationalize these observations. 

2.2 ID vs. 2D blisters 

We first recall the case of the one- dimensional blisters [23J, where a buckled elastic 
line adheres on a substrate. A balance between adhesion (which favors large areas of 
contact) and bending (which is minimized for blisters of large widths) energies dictates 
the shape of the observed blisters. Indeed, within the limit of small amplitudes, the 
adhesion energy density scales as £ 7 ~ 7, while the bending energy density is proportional 

2 Such wrinkles are reminiscent of secondary cracks than can appear due to tensile stresses parallel to 
a main crack at its tip in the presence of a locally weak interface. They act as crack-stoppers, as described 
in the Cook-Gordon mechanism 1221. 
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Figure 2: Typical experimental procedure (/i = 30 /im, AZo = 0.97 mm). (a) The profile of 
the blister is determined by scanning the blister with a tilted laser sheet, while a standard 
image correlation technique is used to estimate the local displacement along the blister: 
the white line on the right corresponds to the intersection of the laser sheet with the 
thin sheet; the white arrows represent the in-plane displacement field, (b) The transverse 
displacement along the blister is found to vary linearly with the distance to the edge X e d ge , 
which allows for a precise definition of the length of the blister L and the position of its 
tip. (c,d) Height and width as a function of the distance from the tip x = L — X e d ge . The 
experimental data are fitted by power functions with respectively 2/3 and 1/3 exponents 
(dashed lines: conical situation). 



to £b ~ B(d/\ 2 ) 2 , where B = Eh? / '12(1 — v 2 ) is the bending modulus of the film. At 
equilibrium, both energies are of same order, which leads to a constant curvature of the 
crest: 

X 2 /d ~ sjBfi = L ec . (1) 

The characteristic length scale L ec compares surface tension (or adhesion) to bending 
stiffness and has been refereed to as elasto- capillary length |24| [25| I26j . 

We compare this result with our 2D experiments conducted with samples of different 
thicknesses h and initial transverse displacements AZo m figure [4] The width A is found 
proportional to d 1 ^ 2 , with a different prefactor depending only on the thickness h (inset). 
Moreover, all the experimental data are found to collapse into a single power law A ~ 
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Figure 3: Profiles of the successive sections of a same blister represented in the normalized 
coordinates (z / cIq) / (x / X) 2 / 3 vs. (y / Xq)/(x/L) 1 ^ 3 (Equations (|4|). Solid line corresponds 
to 1/2[1 + cos(2ttx)]. 



\/dL ec . This result suggests that all the scanned profiles of the blisters behave as if they 
were isolated ID blisters. 



A/4 C 



o.i 





















«< is 






































o J^r"° 






— 


2 






PP 15 urn □ 






PP 30 Mm * 






PP 50 urn o 



0.01 



d/L ec 



o.i 



Figure 4: Normalized width X/L ec as a function of normalized amplitude d/L ec for experi- 
ments conducted with samples of different thicknesses and initial transverse displacements. 
The solid line corresponds to a power law fit with an exponent 1/2: X^/d — 7.3L ec (see 
Equations ([I]) and (|6])). Inset: dimensional data. Colored symbols are used to differentiate 
the different experiments conducted with sheets of the same thickness. 



We have seen that the transverse displacement follows a linear law: 

Al = Al x/L. 



(2) 
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Along a cross section of the blister conservation of length imposes that the total arc length 
of the blister is equal to the sum X(x) + Al(x), which leads (in the limit of small amplitudes) 
to: 

d(x) 2 ~ X{x)Al(x). (3) 

The combination of this last geometrical relation with the expressions for the transverse 
displacement and the curvature finally gives the variations of the amplitude and the width 
of the blister with the distance from the tip: 



d/Ar /J L 



2/3 T 1/3 



(x/L) 



2/3 



and X/All /3 L 2 J C 3 ~ (x/L) 



1/3 



(4) 



which is in agreement with the self-similar profiles displayed in figure [3] The data corre- 
sponding to samples of different thicknesses and initial transverse displacements are super- 
posed in figures^ and [5b. Although scattered, the data collapse into the predicted master 
curves over two orders of magnitude (measurements obtained close to the tip tend to be 
noisy due to the very small amplitude involved and to the presence of secondary wrinkles 
in the vicinity of the tip). 





(a) 



(b) 



Figure 5: (a,b) Rescaled height and width as a function of the dimensionless distance from 
the tip (Equations Q). Solid lines correspond to power law fits with respective exponents 
2/3 and 1/3. 



2.3 Stretching energy in 2D blisters 

We have seen that the shape of a blister can be viewed as a continuous succession of 
profiles corresponding to one-dimensional blisters. However, this shape should correspond 
to the minimum of the global energy of the system. Contrary to the case of an elastic 
line, the bending energy of a plate involves the two principal curvatures. However, if the 
curvature of the crest line is not too high (d"(x) <C d/X 2 ), the bending energy of the blister 
reduces to the ID estimate. 

More importantly, since the thin sheet is curved along the crest of the blister in its two 
principal directions, it does not conserve its initial zero gaussian curvature (the product of 
both principal curvatures). According to Gauss' Theorema Egregium |27j . this change in 
gaussian curvature implies stretching or compression in the thin sheet. Stretching energy is 
usually so dominant in thin objects, that it leads to the formation of crumpling singularities 



to minimize it [32] . One way to estimate the stretching energy is to compare the actual arc 
length of the crest to the length of the non stretched ideal cone (figure |2^) . If the profile 
is given by d(x), the curvilinear length is obtained from the integral f Q L (1 + d 1 (x) 2 ) l / 2 dx, 
while the non stretched length is simply L(l + (d®/ L) 2 ) 1 / 2 . If d(x) follows a power law, the 
typical strain is thus given by e ~ (do/L) 2 . We finally obtain an estimate for the stretching 
energy density inside the blister: £ s ~ Eh(do/L) 4 . We expect the description of the blister 
as a succession of ID profiles to be valid in the limit where £ s is negligible when compared 
to £y (or £(,), which corresponds to: 

A«(£) 1/4 i. (5) 

In our experiments both (j/Eh) 1 ^ and do/L are typically on the order of 2%. Nevertheless 
all numerical prefactors have been omitted in our description, and the stretching energy 
may have been largely overestimated. In order to check the validity of these estimates, nu- 
merical simulations were carried out with the Surface Evolver software |28j . The physical 
parameters selected for the simulation correspond to a typical experiment: E = 2000 MPa, 
i/=0.4, h = 50 /um, A/o = 1 mm and 7 = 22.4mN.m~ 1 (corresponding to (-f/Eh) 1 ^ ~ 2%). 
The plate is discretized into finite elements with an elastic energy corresponding to their 
deformation. Adhesion energy is represented by a potential decaying exponentially with 
the distance from the rigid plate (we checked that the characteristic decay length of the 
potential does not influence the final result, as long as it remains small). Surface Evolver 
minimizes the energy of the system using alternately gradient and conjugate gradient algo- 
rithms. Stochasticity is introduced by "jiggling" the position of vertices to escape from local 
minima of energy. Mappings of the bending and stretching energy densities are displayed 
in figures |6^ and [6)3. The maximal values in the blister, respectively 4.10 -5 mJ.mm -2 
and 3.10 _7 mJ .mm for bending and stretching energy densities, are consistent with the 
scalings B/L 2 C ~ 10~ 5 mJ.mm" 2 and Eh(d$/ L) ~ 10~ 6 mJ.mm -2 . These results confirm 
that stretching energy along the blister is indeed negligible in this regime. We also checked 
that the numerical profiles obeyed equations Q. 

In the opposite case, where (p| is not satisfied, stretching energy plays a dominant 
role, and the previous scalings (HF are no longer valid. As a consequence, the shape of 
the blister that minimizes the energy of the system is a developable cone. This limit is 
discussed in details by Chopin et al. [21j. Indeed, a typical value for do in their liquid blis- 
ter test is on the order of 5 mm, while the product L^/Eh) 1 ^ ranges between 1 and 5 mm. 

From a quantitative point of view, the prefactor in equation ([I]) can be found by using 
an exact expression for the the elastic energy of the blister, or equivalently by using the 
fact that the curvature at the adhesion point (and at the crest) is equal to 2/L ec [26J. 
Equation ([I]) becomes finally |14] : 

X 2 /d = ir 2 L ec . (6) 

However, fitting the experimental data (figure Q results in a slightly different prefactor 
X 2 /d = 7.3L ec . This mismatch is not due to an effect of the stretching strain. Indeed 
the gaussian curvature of the crest of the blisters appears larger than in the case of a 
collection of independant one-dimensional blisters, which leads to an higher stretching 
energy. The origin of the discrepancy certainly lies in the experimental preparation of the 
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Figure 6: Numerical simulation using the Surface Evolver software: E = 2000 MPa, 
v = 0.4, h = 50 /jm, AIq = 1mm and 7 = 22.4mN.m~ 1 . Maps of bending (a) and 
stretching (b) energy densities confirm that the stretching energy inside the blister is neg- 
ligible (dashed line: blister's edge). 



blister. During the evaporation of the wetting liquid layer, temporary large depression 
leads to smaller, more curved blisters. This might lead to some plastic deformation in 
the material (plasticity is expected in polypropylene for a radius of curvature smaller than 
lOO/i). For the 15 fi m sheets, the curvature expected in the blister is 3mm, which is not 
very far from lOO/i = 1.5 mm. However no clear signs of plasticity have been observed in 
the experiments. 

2.4 Thin sheets on soft substrates under compression 

In this paper we focus on the situation of an adhesive sheet deposited on a rigid passive 
substrate, and blisters appear because of imposed lateral displacements (controlled, as 
in section 2, or randomly in section 3). But another common situation is found when 
thin films deposited on a substrate delaminate under compressive stresses induced by 
the deposition process, thermal expansion, chemical swelling, or mechanical compression 
of the substrate. For instance, a recent study motivated by potential applications to 
stretchable electronics focused on the formation of one-dimensional blisters observed when a 
soft substrate covered with a thin film is uni-axially compressed |14j . Prior to delamination, 
wrinkles are commonly observed as the result of the balance between the bending energy 
of the film and the stretching energy of the substrate |5] . The wavelength of such wrinkles 
scales as A ~ h(E / Eg) 1 / 3 , where E and E$ are the Young moduli of the film and the 
substrate, respectively. However we focus here on delamination where the film detaches 
from the substrate. In this case, three physical ingredients dictate the buckling behavior 
of the thin layer: adhesion, bending stiffness and, in addition, the compression of the 
substrate. Although the condition to delaminate relies on a combination of these three 
ingredients, this study shows that the shape of one-dimensional blisters is set by the balance 
between adhesion and bending rigidity, identical to the case of a thin sheet able to slide 
on a rigid plate (the transverse displacement depends however on the Young's modulus of 
the substrate). 
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We assume here that the possible formation of wrinkles has a negligible impact while 
the width of the blisters is large in comparison with the wavelength of these wrinkles. 
We thefore expect that our results on angular blisters (eqn. [T]) extend to the case of the 
delamination of thin films from soft substrates if condition eq. [5] is satisfied. To validate 
this assumption, a sheet of polypropylene (E = 2.1 ± 0.2 GPa, h = 30 ^m) is uniformly 
laid on a block of thickness hs = 5 cm made with vinylpolysiloxane (Zhermack) of Young's 
modulus Es — 220 kPa. The thin film adheres to the substrate through Van der Waals 
interactions. The block is obliquely compressed with non-parallel clamps, which induces 
the propagation of one or multiple angular blisters in the direction perpendicular to the 
compressive load (sketch in the insert of figure [7J . The shape of a typical blister reported 
in figure [7] follows the same form A ~ l}J^d l l 2 as in the rigid substrate configuration. The 
length of the blister L ~ 80 mm is sufficient to fulfill the condition (5). The deduced value 
of the elasto- capillary length L ec = 9.6 mm is in good agreement with the estimation from 
material properties and with the prefactor for the law A ~ L\^d l l 2 found in Vella et. 
al.gQ. 




d {mm) 

Figure 7: Compression of a soft substrate covered with a thin film (insert: sketch of the 
experiment): as in the rigid substrate configuration, the curvature of the crest is found 
constant (A ~ Le^d 1 / 2 ) as prescribed by Equation [gJ 

In conclusion to the simplified case of a single angular blister with a small slope, the 
radius of curvature of the crest is found constant and proportional to L ec . Measuring this 
radius of curvature thus allows to estimate the elasto -capillary length, a typical length scale 
which compares bending stiffness and adhesion, (or one of these quantities if the other one 
is known). The robustness of this result to simultaneous random blisters is probed in the 
following section. 

3 Multiple random blisters 

We present in this section the generalized case of multiple blisters obtained as a thin 
sheet is randomly laid on an adhesive substrate, and also discuss the shape of the tips of 
the blisters. 

Contrary to the model experiment with a single angular blister, the random blisters 
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may be closed and trap a certain volume of air between the thin sheet and the substrate. 
We assume here that the pressure is small enough and do not influence the shape of the 
blister 4 {] In the case of high pressure, the contact line would tend to be circular |!6| 117] . 
in contrast with the elongated shapes observed in our experiments. 

3.1 Curvature profiles 

A rigid plate is preliminary coated with a thin layer of vinylpolysiloxane, which pro- 
vides adhesion through van der Waals interactions. The elasto- capillary length is pre- 
liminary measured from a one-dimensional blister and gives L ec ~ 22 mm. Multiple blis- 
ters are randomly formed afterward by depositing carelessly a thin sheet of polypropy- 
lene (E = 2.1 ± 0.2 GPa, h = 30 /mi) over the plate. The height profile of the whole 
sheet is determined with a versatile optical profilometry technique based on the deforma- 
tion of a fringe pattern projected on the surface |29[ I30j . The field of mean curvature 
C m = 1/2(1/ R\ + I/R2), where R\ and R2 are the principal radii of curvature, is readily 
deduced from the topography with a standard MatLab routine. Although blisters with 
different heights and slopes are observed (figure [8k) , the curvature of the crests remains 
almost uniform for the different blisters (figure pp) . In the particular case illustrated in 
figure [8| we find C m ~ 0.04 mm -1 , which corresponds to L ec = 1/C m ~ 25 mm in good 
agreement with a preliminary measurement with a one-dimensional blister. 

One blister in the middle of the figure |8j3 seems however significantly different from 
the others: its crest is not straight but slightly curved, with a typical radius of curvature 
p sketched in figure [8};. As friction prevents re-arrangements of the adhering patches, such 
curved blisters tend to appear randomly depending on the way the film has been deposited. 
This curvature induces gaussian curvature, and therefore an additional source of stretching 
energy. An estimate of this gaussian curvature K is given by the product of the typical 
curvature of a section of the blister (i.e. 1/L ec , as shown in the previous section) by the 
curvature of its axis projected on the normal to the sheet (inclined with a typical angle 
d/X) |27| : K ~ (l/L ec )(d/\p). To compute the typical strain e, we use the fact that 
Gauss' Theorema Egregium can be written as K ~ d 2 e |31| I32j . This scaling can also be 
deduced from the second Foppl-von Karman compatibility equation A 2 + EK = [33J, 
where <f> is the Airy stress function (with d 2 (f) ~ Ee) and E the Young modulus of the 
material. Since the width A corresponds to the typical length of variation of the strain in 
the blister, the strain can be estimated as e ~ Xd/pL ec . Note that in the case of straight 
angular blisters, the same derivation leads to K ~ d 2 /X 2 x 2 , and thus to an estimate of 
the strain e ~ d 2 /x 2 , which is consistent with the equation pM The density of stretching 
energy of a curved blister is thus on the order of EhX 2 d 2 / p 2 ]^ 2 ^ Neglecting the effect of 
the curvature of the blister is only valid if EhX 2 d 2 / p 2 L 2 c < 7, ie. /) > Xd/h. For the 
experiment shown in figure [HJ the critical radius beyond which this additional stretching 
energy can be neglected is on the order of 700 mm, which is higher than the typical radius 
observed for this particular blister (p ~ 50 mm). The additional stretching energy should 
explain the peculiar profile of this curved blister. 

Note finally that in the model experiment described in the previous section, the move- 
ment of the sides of the single blister is a rotation, whose center is located at the tip of 
the blister. In the case of multiple random blisters, each tip of a blister is a local center of 

4 The effect of a pressure p is negligible if the associated energy per unit surface pd is much smaller 
than 7, or equivalently p <C 7/d, where d is the heigh of the blister. 
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Figure 8: (a) Height profile of random blisters formed when a thin sheet is carelessly 
deposited over an adhesive plate, (b) Field of absolute mean curvature \C m \ = \\/2{l/ R\ + 
1/.R2)|, where R\ and R2 are the principal radii of curvature, (c) Sketch of the different 
radii of curvature. 

rotation of the adhered parts of the thin sheet. 

It would be tempting to extrapolate our macroscopic results down to the molecular 
scale of graphene sheets. For instance, Xu et al. [3] recently reported blisters on mono- 
layer graphene sheets deposited on SiC>2 substrates of widths, heights and lengths on the 
order of 10 nm, 3nm and 1 ^m, respectively. We can thus estimate from the ratio X 2 /d 
the corresponding elasto- capillary length L ec ~ 3nm. Assuming a bending stiffness for a 
monolayer on the order of 0.2nN.nm derived from phonon spectrum of graphite, poten- 
tials or ab initio calculations |34l l35| 136] . we find an adhesion energy 5 A7 on the order of 
50mJ.m _2 . [] Although ab initio calculations predict an adhesion energy one order of mag- 
nitude higher |37j . large mismatches have been reported between calculated and measured 
values in similar situations [4j. Moreover, the height of these blisters are on the order of 
the typical decay length of Van der Waals interactions [37] and such interactions cannot be 
neglected to compute the shape of the buckled part of the graphene sheet. Progresses have 
still to be made to produce a reliable estimate of the adhesion energy at those microscopic 

5 When adhesion is provided by a wetting liquid, the work of adhesion is equal to 27. In the case of 
conventional adhesion, the adhesion energy A7 is thus equivalent to 27 in the definition of L ec . 
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scales, but the observation of blisters may constitute a valuable method. 



3.2 Radius of the tip 

We finally focus on another feature of these small amplitudes blisters: the shape of the 
tip. The scalings previously shown for the width and amplitude of blisters (equation [4]) are 
no longer valid near the tip, where they would lead to a diverging curvature of the crest line. 
Indeed, since d(x)/x is proportional to x -1 / 3 , the local version of condition ^ is always 
violated close to the tip of the blister. As a matter of fact, any smooth regularization of the 
tip necessarily involves non-zero gaussian curvature, producing in-plane deformations. We 
thus assume that all three energies (bending, stretching, and adhesion) are of same order 
in this region. If the tip ends with a circular contact line (radius R), the stretching energy 
density involved in the spherical cap of amplitude d is proportional to Eh(d/R) 4 following 
an estimate similar to equation [HJ while the bending energy density scales as Eh 3 d 2 / R 4 . 
Imposing that both energies are comparable to the adhesion energy 7 leads to d ~ h and 
R 2 /d ~ L ec . We finally find that the radius of the tip is given by: 

R ~ ^hZ7 c . (7) 

A more quantitative description of the shape of the tip would require solving Foppl-von 
Karman equations numerically. We nevertheless propose here a simplified model where we 
assume that the shape of the tip is the half of an axisymetric cap of radius R and height 
d. The Lagrangian |38j (or more commonly, the energy) of the cap can be approximated 
I'.v'O 

£=— ^+ R2 +TrRh + P(d-d Up ), (8) 

where a is a numerical constant (a ~ 0.07) and /3 the Lagrange multiplier that imposes 
the continuity of height (at d = du p ) between the tip and the remaining part of the blister. 
Minimizing C with respect to R, d and ft yields: 

R = . (9) 

The amplitude du p has to be related with the height beyond which an axisymetric blister 
is no longer stable, so, as pointed out in [21|, du p ~ Sh, where 5 is a numerical prefactor: 

R = {AS 2 + 12a(l - v 2 )b A f^^fhL^ c . (10) 



We find here the prefactor for the proposed scaling law (equation (|7j)), but it depends on 
the additional parameter 5. Experimental data obtained using films of various polymeric 
materials show that the radius of the tip varies linearly with \/hL ec and can be approx- 
imated by R = 3.2y/hL ec (figure [9]). Despite the variations on both the measurements 
of R and the elasto- capillary length 7 , as shown by the error bars in figure [9J this result 
appears significant knowing that the thicknesses used range from 15 /im to 1 mm. [] The 

6 See paragraph 100 in [38] . valid for deflections of order h. 

7 The value for adhesion energy of the polymer on the vinylpolysiloxane depends if the measurement 
corresponds to sticking a lamella of polymer on the substrate (lower adhesion), or to debounding (higher 
adhesion). The horizontal error bars in figure[9]have been determined with these extremes values, obtained 
by measuring the shape of one-dimensional blisters. 
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prefactor found experimentally yields an estimate of the parameter 5 ~ 3.1. Despite the 
fact Chopin et al. noticed that the threeshold for the instability of the axisymetric blister 
was not reproducible in their experimental set-up, we can estimate from their data a crit- 
ical amplitude d c ~ 50h. However, in the liquid blister case, the shape loses its symmetry 
and becomes triangular after the onset of the instability. Further work is still necessary to 
characterize completely the shape of the tips of these delamination blisters. 

We conclude that while the ratio X 2 /d provides an estimate of L ec , measuring R gives an 
additional information on the thickness of the sheet. Simple topographic characterization 
of random blisters finally allow to access both relevant length scales of the problem. 

8 
7 
6 
5 

S 

S 4 
^ 3 
2 
1 


0.5 1 1.5 2 2.5 

V hL„. {mm) 

Figure 9: Radius of the tip R as a function of \JhL ec for sheets of various polymeric mate- 
rials and thicknesses ranging from 15 fim to 1 mm (Equation (J7|). R can be approximated 
by R = 3.2y/hL^ c . 
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4 Conclusion 

Although the formation of blisters on thin sheets deposited on a substrate is usually 
viewed as an undesirable defect, a simple geometrical characterization of low amplitude 
blisters provides estimates of both the thickness of the sheet and an elasto- capillary length 
scale L ec = \jBj^ that compares bending stiffness B and adhesion 7. The radius of 
curvature of the crest is found uniform through the whole sample and proportional to 
L ec while the radius of the tip is proportional to ^hL ec . Scaling down our macroscopic 
experiments may allow to estimate the mechanical properties of thin polymeric films or 
ultimately of molecular films such as graphene sheets that are difficult to assess. 
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